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Heavy Flavour Physics at CDF 

Gavril Giurgiu on behalf of the CDF collaboration 

Phyiscs and Astronomy Department, The Johns Hopkins University 

3400 N. Charles Street, Baltimore, MD 21218, U.S.A. 

The CDF detector at Fermilab has accumulated more that 3 fb _1 of data which enables unprecedented studies 
of heavy flavor hadron properties. We present recent CDF measurements of mass and lifetime of the B c meson 
as well as lifetime, mixing and CP violation properties of B s mesons. 



1. Introduction 

We review recent measurements of B c and B s 
meson properties performed with the CDF de- 
tector in y/s = 1.96 TeV pp collisions at the 
Tevatron. Due to the excellent performance of 
the Tevatron the CDF experiment has recorded 
3.5 fb^ 1 of data by June 2008. The measure- 
ments presented here use between 1 fb~ x and 
2.4 fb^ 1 . The CDF detector is described in detail 
in Ref. [J. 

At the Tevatron, B hadrons are mostly pro- 
duced in pairs. The main bb production mecha- 
nism is flavor creation through gluon fusion. The 
bb production cross section of ~ 30 [ib [I] is 
large compared to B-factories which enables rich 
B Physics programs at the Tevatron. However, 
the Tevatron bb cross section is orders of magni- 
tude smaller than the total inelastic cross section 
of ~ 50 mb. For this reason, CDF employs trig- 
gers that select events with signatures specific to 
various B decays. The measurements presented 
here uses data selected by the di-muon trigger, 
geared towards B — > J/ipX decays and by the 
displaced track trigger used for measurement of 
B s lifetime in B s —> D s irX decays. 

2. Mass and Lifetime of B c Mesons 

The B c meson is unique as it contains two 
heavy quarks: b (bottom) and c (anti-charm). 
Different theoretical models predict the B c mass 
around 6.3 GeV/c 2 . In particular, non-relativistic 
potential models [2] predict the range 6247 — 
6286 MeV/c 2 . Lattice QCD models [3] predict 



6304 ± 12±q 8 MeV/c 2 , while similar results are 
determined in perturbative QCD calculations [4]. 

The lifetime of the B c meson is expected to be 
~ 0.5 ps [5] which is about a third of the typi- 
cal light B meson lifetime of ~ 1.5 ps. The short 
lifetime is explained by the fact that either the 
b or the c quarks can decay weakly through a 
W boson. When the b quark decays, the final 
state contains a charmonium particles while the 
decay of the c quark results in a B meson final 
state. Additionally, the b and c quarks can an- 
nihilate, for example into a lepton and neutrino 
final state. The measurement presented here use 
B c — > J/ipli/X decays, where the lepton is either 
a muon or electron. 

2.1. Measurement of B c Mass 

The mass of the B c meson was measured 
for the first time by the CDF experiment in 
Run I [B] using inclusive semilcptonic decays 
B c — > J/i/jIvX. Since the B c meson was 
not fully reconstructed due to the missing neu- 
trino, the result had large uncertainties 6400 ± 
390(stot.) ± 130(sysi.) MeV/c 2 . 

The most precise measurement of the B c mass 
was recently performed by the CDF experiment 
using 2.4 fb^ 1 of data [7] using the fully recon- 
structed decay mode B c — > J/ijj tt. The advan- 
tage of using fully reconstructed decays is that 
the mass of the decaying particle can be measured 
precisely by fitting the invariant mass distribu- 
tion. The measurement was performed using data 
collected by the CDF J/ip trigger. This trigger 
requires two muons of opposite charge, transverse 
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momentum larger than 1.4 GeV/c and invariant 
mass close to the J/tp mass. The data sample 
used in this analysis contains ~ 17 million J/tp 
candidates reconstructed with an average mass 
resolution of 13 MeV/c 2 . 

Using this initial J/tp sample, both B + — ► 
J/tp K + and B c — > J/tp tt + decays are recon- 
structed. The B + sample is used to determine 
an optimal set of selection requirements that are 
subsequently used for the B c selection. This pro- 
cedure accounts for the expected shorter B c life- 
time by using only B + mesons with decay time in 
rages typical to expected B c lifetime: 80 < ct < 
300 fim. 

The invariant mass distribution of the B c can- 
didates is shown in Figure [1] The mass of the 
B c meson is determined using an un-binned log 
likelihood method. The signal probability den- 
sity function (PDF) is described by a Gaussian 
function with a width proportional to the esti- 
mated mass uncertainty estimated from track un- 
certainties. The background PDF includes both a 
combinatorial component as well as a component 
from B c — ► J/tp K + Cabibbo suppressed decays. 
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Figure 1. B c invariant mass distribution recon- 
structed from B c — > J/tp 7r decays. 



The total B c signal yield is estimated to 108 ± 
15 candidates with a statistical significance of 
8 standard deviations. The measured mass is 
6275.6 ± 2.9 MeV/c 2 . The largest systematic 
uncertainty of 2.2 MeV/c 2 is due to the knowl- 



edge of the scale factor on mass uncertainties. 
The total systematic uncertainty which includes 
effects due to misalignment, momentum scale 
and Cabibbo suppressed B c — > J/tp K + decays 
is estimated to 2.5 MeV/c 2 , comparable but 
slightly lower than the statistical error, leading 
to the most precise measurement of the B c mass: 
6275.6 ± 2.9{stat.) ± 2.5(syst.) MeV/c 2 , in good 
agreement with the corresponding DO measure- 
ment [§] 6300 ± U{stat.) ± 5(syst.) MeV/c 2 

2.2. Measurement of B c Lifetime 

The lifetime of the B c meson is measured by the 
CDF experiment using partially reconstructed 
semileptonic B c — ► J/tp (/i/i) I v X decays [TU] in 
1.0 fb^ 1 of data. The lepton I can be either a 
muon or an electron. The data sample is se- 
lected by the same J/tp trigger used for the B c 
mass measurement. The lifetime analysis starts 
with a sample of ~ 5.5 million J/tp — > /i + \i~ 
candidates. A third lepton is then required to 
come from the same vertex as J/tp. The main 
challenges in this measurement are the partially 
reconstructed momentum of the B c meson and 
multiple backgrounds. 

Since the momentum of the B c meson is not 
fully reconstructed, the observed pseudo proper 
decay time ct* = " l ^ v }f£t) ^ ^ s statistically cor- 
rected by a factor K — Exi^M*} determined from 
simulation. Here, m is the mass of the B c me- 
son [5], pr(J/tpl) is the reconstructed transverse 
momentum of the J/tp and lepton, px(B c ) is the 
true transverse momentum of the B c meson deter- 
mined from simulation and L xy is the transverse 
B c decay length projected on the reconstructed 
transverse momentum. 

In this measurement the mass of the B c meson 
is not reconstructed due to the missing neutrino. 
Instead, the J/tp I mass is used. However, this 
quantity has a wide distribution and sidebands 
cannot be used to estimate backgrounds as in 
fully reconstructed decays. Different background 
contributions are determined using data-driven 
methods when possible and simulation otherwise. 
The main sources of backgrounds are: fake lep- 
tons, fake J/tp, "66 sources", e + e~ from either 
photon conversions or decays of ir° and r\ mesons 
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and prompt J/ip events. The fake lepton back- 
ground is produced by a real J/ip and a ran- 
dom track which fakes a muon or an electron 
and forms a J/ip I vertex that satisfies the vertex 
quality requirements. The fake lepton fractions 
are determined using real pions and kaons from 
D*+ -> D° 7T+ with D° -> 7T+ K- and real pro- 
tons from A — > p + ir~ , where the real pions, kaons 
and protons are misidentified as leptons. The fake 
J/ip background is modeled using the J/ip mass 
sidebands, bb backgrounds are events with a real 
J I ip from one b quark and a real lepton from the 
other b anti-quark in the event. This background 
is modeled using Pythia simulation where the 
fractions of the three main production mecha- 
nisms, flavor creation, flavor excitation and gluon 
splitting are tuned to data using the angular sep- 
aration between the J/ip and the lepton as dis- 
criminating variable. e + e~ backgrounds, specific 
to J/ip e mode only, are suppressed by trying to 
identify both the positron and electron from pho- 
ton conversion or light meson decay. The residual 
e + e~ pairs are modeled using Pythia simula- 
tion. The background fraction from prompt J/ip 
is modeled by a Gaussian function in the pseudo 
proper decay space and determined directly from 
the un-binned likelihood fit. 

The B c lifetime is measured separately in J/ip e 
and J/ip fi channels using un-binned likelihood 
methods. The fit projections of the pseudo proper 
decay time distributions in the J/ip fi channel is 
shown in Figure [2j 

The corresponding lifetimes are: cr M = 
n^.ltf 7 %{stat.) and cr e = 121. 7 tlt' 3 {stat.) fim 
The main sources of systematic uncertainties 
come from our understanding of the decay time 
resolution function (3.8 fim) and from relative 
fractions of bb production mechanisms in Pythia 
simulation. The total systematic uncertainty 
is 5.5 fim, leading to a combined measurement 
of the B c lifetime of ct m = 142.5i^'j (stat.) ± 
5.5(syst.) fim. This agrees well with a similar 
recent measurement performed by the DO ex- 
periment 12] in the B c — ► J/ipfi v X channel 
ct= 134A±ll^(stat)±9.6(syst). Figure 1 shows 
a summary of the experimental results of the B c 
lifetime. CDF and DO measurements agree well 
with the theoretical predictions and with each 
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Figure 2. Pseudo proper decay time fit projec- 
tions in B c — > J/ip fiv X channel. 



other and have similar uncertainties. 

3. Neutral B s System 

A B° B meson is a bound state composed of an 
anti-bottom quark b and a strange s quark. The 
time evolution of a mixture of the B® and its an- 
tiparticle B°, a(t)\B°) +b(t)\B°), is given by the 
Schrodinger equation 

<!(;)-("-'§)(:)■ 

where M and T are 2x2 mass and decay matri- 
ces. The mass eigenstates B\ and B® H are liner 
combination of the flavor eigenstates and B® 
and are obtained by diagonalizing the mass op- 
erator. The mass difference Am s is proportional 
to the B s mixing frequency recently measured by 
both CDF [J5] and DO [T3] experiments and found 
to be in good agreement with Standard Mmodel 
(SM) predictions |17) . Moreover, the mass eigen- 
states have different decay widths (lifetimes) Tl 
and Th- The average decay width is defined as 
T = (r L + T H )/2 and the decay width differ- 
ence is defined as Ar = Tl — Th- The decay 
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Figure 3. Summary of the B c lifetime experi- 
mental results. 



width difference Ar = 2|Ti2 1 cos(0 s ) is sensitive 
to new physics (NP) effects [15116117] that affect 
the phase <p s = arg(— M 12 /r 12 ), where T 12 and 
M12 are the off-diagonal elements of the mass 
and decay matrices. New physics will increase <p a , 
so Ar would be smaller than the SM prediction. 
Since the SM phase </>f M is predicted to be very 
small (~ 0.004) [IT], in a new physics scenario 
with large contribution to (f> s one could approx- 
imate (j) s = <j) S s M + (f>^ p w p . If real, this 
new physics phase is accessible in B s — > J/ip <p 
decays. In these decays one can measure the 
CP violation phase f3 s which is the relative phase 
between the direct decay amplitude and mixing 
followed by decay amplitude. In SM this phase 
is defined as /3 S W = arg(-^$|) QU, where 

Vij are the elements of the CKM quark mixing 
matrix. Global fits of experimental data tightly 
constrain the CP violation phase to small val- 
ues in the context of SM /3f M w 0.02 [18]. The 
presence of new physics could modify this phase 
by the same quantity <p^ p that affects the <pf M 
phase. New physics could contribute significantly 
to the observed p3 s phase [15116117) expressed as 



2(3 8 = 2f3g — <p s . Assuming that new physics 
effects dominate over the SM phase, we can ap- 
proximate 2{3 S « ~(t>^ p ~ —<p s - 



3.1. B s Lifetime and Width Difference 

A sample of ~ 2500 signal B s — ► J/ip cf> events 
collected with the CDF J/ip trigger in 1.7 fb^ 1 of 
data was used to measure the B s lifetime t s and 
decay with difference Ar s [19] . The B s invariant 
mass distribution is seen in Figure [4] which shows 
good signal to background ratio. The B s meson 
is a spin particle which decays in two spin 1 
particles J/ip and (p. The total angular momen- 
tum in the final state can be 0, 1 or 2. The states 
with angular momentum and 2 are CP even, 
while the state with angular momentum 1 is CP 
odd. One can statistically separate the CP even 
and CP odd states using their different angular 
distributions. In the case of CP conservation the 
mass eigenstates are also CP eigenstates. It is 
interesting to note that since in SM the CP vio- 
lation phase in the B s system is very small, the 
mass eigenstates are CP eigenstates with good 
approximation. 

The B s lifetime and decay width difference are 
extracted using an un-binned likelihood method 
in the space of mass, decay time and angles of the 
final state particles. There are three angles p = 
{cos 9t,<Pt, cos ipr} that completely define the di- 
rections of the four particles [i + , yT , K + , K~ in 
the final state. These angles are defined in the 
transversity basis introduced in Ref. [20]. Fig- 
ure [5] shows the B s decay time distribution to- 
gether with the fit projection. The largest sys- 
tematic uncertainties on lifetime are due to de- 
cay time resolution model and alignment of the 
silicon detector. The largest systematic on Ar 
is due to about 3% of B° — > J/ip K*° decays re- 
constructed as B s — * J/tp (j>. The final results 
are: ct = 1.52 ± 0.04(staf.) ± 0.02(syst.) ps and 
Ar = 0.076t° ^ll(stat.) ± 0.006(syst.)p S - 1 . The 
fractions of each of the three final state ampli- 
tudes are also measured and together with the 
lifetime and width difference are found to be 
in agreement with the similar measurement per- 
formed by the DO experiment [2T] , 
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Figure 4. B s invariant mass from B s — > J/ip < 
decays. Fit projections are superimposed. 



Figure 5. B s decay time from B s — > J/tp(j) de- 
cays. Fit projections are superimposed. 



3.2. CP Violation in B s System 

The CP violation phase j3 s in B s — ► J ftp cj> de- 
cays was measured for the first time by the CDF 
experiment using 1.4 fb^ 1 of data [55]. The anal- 
ysis uses a sample of ~ 2000 signal events. To en- 
hance our sensitivity to the CP violation phase, 
we identify the flavor of the B s or B s meson 
at production by means of flavor tagging. Two 
independent types of flavor tags are used, each 
exploiting specific features of the production of 
b quarks at the Tevatron. The first type of fla- 
vor tag infers the production flavor of the B s or 
B s meson from the decay products of the b hadron 
produced by the other b quark in the event. This 
is known as an opposite-side flavor tag (OST). 
The second type of flavor tag identifies the flavor 
of the reconstructed B s or B s meson at produc- 
tion by correlating it with the charge of an asso- 
ciated kaon arising from fragmentation processes, 
referred to as a same-side kaon tag (SSKT). The 
average dilution D, defined via the correct tag 
probability P = (1 + D)/2, is (11 ± 2)% for the 
OST and (27 ± 4)% for the SSKT. The mea- 



sured efficiencies for a candidate to be tagged are 
(96 ± 1)% for the OST and (50 ± 1)% for the 
SSKT. 

An unbinned maximum likelihood fit is per- 
formed to extract the parameters of interest, 2(3 a 
and Ar, plus additional parameters referred to as 
"nuisance parameters" which include the signal 
fraction f s , the mean B s width r = (Tl +Ph-)/2, 
the mixing frequency Am s , the magnitudes of the 
polarization amplitudes \A \ 2 , \A\\\ 2 , and |Aj_| 2 , 
and the strong phases 6\\ = arg(Aj|Ao) and 5± = 
aig(A* L Ao). The fit uses information on the re- 
constructed B s candidate mass m and its uncer- 
tainty cr m , the B s candidate proper decay time 
t and its uncertainty cr*, the transversity angles 
p and tag information T> and £, where T> is the 
event-specific dilution and £ = { — 1, 0, +1} is the 
tag decision, in which +1 corresponds to a can- 
didate tagged as B s , —1 to a B s , and to an 
untagged candidate. 

The time and angular dependence of the signal 
PDF P s (t, p, £, \T>, at) can be written in terms of 
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two PDFs, P for B s and P for B s , as 



P s (t,p^\V,a t ) 



2 

1-SV 



P(t,p\at)e(p) 
P(t,p\a t )e(p), (2) 



The detector acceptance effects on the transver- 
sity angle distributions e(p) are modeled with 
B s — ► J/^i simulated data. The time and angu- 
lar probabilities for B s can be expressed as 

P(t,p) cx \A Q \ 2 T + h(p) + \A^T + f 2 {p) 

+ |A ± | 2 r_/ 3 (p) + |A ll || J 4 ± |w + / 4 (p) 

+ |Ao||A|||co S (d||)T + / 5 (p) 

+ \A \\A x \V+f 6 (p), (3) 

where the functions f\{p) ■ ■ ■ fe(p) are defined in 
Ref. 19 1 . The probability P for B s is obtained 
by substituting U + — > U- and V + — > V_. The 
time-dependent term 7± is defined as 



-rt 



[cosh(Art/2) =F cos(2/3 s ) sinh(AR/2) 
=F ?ysin(2/3 s )sin(Am s t)] , 



where 77 — +1 for P and —1 for P. The other 
time-dependent terms are defined as 



U± = ±e 



-rt 



V± = ±e 



-rt 



x [sin(<5j_ — 5\\) cos(Am s t) 

— cos(S± — Sn ) cos(2/3 s ) sin(Am s i) 

± cos(<S ± - 6\\ ) sin(2/3 s ) sinh(Art/2)] 

x [sin(^x) cos(Am s t) 

— cos(S±) cos(2/3 s ) sin(Am s t) 

± cos(S_l) sin(2/3 s ) sinh(An/2)] . 



The time-dependence is convolved with a Gaus- 
sian proper time resolution function with stan- 
dard deviation cr t , which is adjusted by an overall 
calibration factor determined from the fit using 
promptly decaying background candidates. The 
average of the resolution function is 0.09 ps, with 
a root-mean-square deviation of 0.04 ps. 

Possible asymmetries between the tagging rate 
and dilution of B s and B s mesons have been stud- 
ied with control samples and found to be statis- 
tically insignificant. We allow important sources 
of systematic uncertainty, such as the determi- 
nation of overall calibration factors associated 
with the proper decay time resolution and the 



dilutions, to float in the fit. The mixing fre- 
quency Am s — 17.77±0.12 ps -1 is constrained in 
the fit within the experimental uncertainties 13J . 
Systematic uncertainties coming from alignment, 
detector sculpting, background angular distribu- 
tions, decays from other B mesons, the modeling 
of signal and background are found to have a neg- 
ligible effect on the determination of both Ar and 
j3 s relative to statistical uncertainties. 

An exact symmetry is present in the signal 
probability distribution, as can be seen in Eq. ([3]), 
which is invariant under the simultaneous trans- 
formation (2(3 S -> 7T - 2/3 s , AT -> -Ar, S\\ -> 
2ir — St\ , and S± — > tt — 5± ) . This causes the likeli- 
hood function to have two minima. Since the log- 
likelihood function is non-parabolic and the two 
minima are barely separated at one standard de- 
viation level, we cannot meaningfully quote point 
estimates. Instead we choose to construct a con- 
fidence region in the 2(3 S — AT plane. 

We use the Feldman-Cousins likelihood ratio 
ordering to determine confidence levels (CL) in 
2(3 S — AT space. The other parameters in the fit 
are treated as nuisance parameters {e.g. B s mean 
width, transversity amplitudes, strong phases). 
To ensure that the obtained confidence regions 
provide the quoted coverage against deviations of 
the nuisance parameters from their values mea- 
sured in our fit to data, we perform pseudo- 
experiments by randomly sampling the nuisance 
parameter space within ±5cr of the fit values and 
confirm coverage of the 68% and 95% confidence 
regions shown in Fig. [6] The solution centered 
inO < 2f3 s < tt/2 and Ar > corresponds 
to cos(<5j_) < and cos(S± — 5\\) > 0, while 
the opposite is true for the solution centered in 
tt/2 < f] s < tt and Ar < 0. Assuming the 
standard model predicted values of 2/3 s = 0.04 
and Ar = 0.096 ps^ 1 [17 , the probability to ob- 
serve a likelihood ratio equal to or higher than 
what is observed in data is 15%. Additionally, we 
present a Feldman-Cousins confidence interval of 
2(3 s , where Ar is treated as a nuisance parame- 
ter, and find that 2/3 s e [0.32, 2.82] at the 68% 
confidence level. The CP phase 2f3 s , AT, T, and 
the linear polarization amplitudes are consistent 
with those measured in Ref. [TH] and in Ref. [3T] . 
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Figure 6. Confidence region in the 2(3 S — AT 
plane, where the standard model favored point is 
shown with error bars [T7] . 



3.3. B s Lifetime in Flavor Specific Decays 

The lifetime of the B s meson is measured by the 
CDF collaboration in flavor specific B s — ► D s it X 
decays [33] using 1.3 of data. The flavor of 
the B s meson at decay is given by the charges of 
the final state particles, hence the name "flavor 
specific" . These events are collected by a trigger 
which requires two tracks with impact parameter 
between 120 /im and 1 mm, specific to hadronic 
B decays. The sample contains ~ 1100 fully re- 
constructed B s — > D s (<j>n) n and ~ 2200 partially 
reconstructed decays. Figure [7] shows the B s in- 
variant mass distribution and details of all con- 
tributions to partially reconstructed decays. 

There are two main challenges in this anal- 
ysis. First, the lifetime distribution of the B s 
meson is biased due to the trigger that requires 
tracks with large impact parameter. This is cor- 
rected by using an efficiency function determined 
from simulation. This method is validated on 
control samples of B° -> D~iz + , B° -> D*-ir+ 
and B + — ► D~n + decays where good agreement 
with the B° and B + world average lifetimes is 



Figure 7. Invariant B s mass distribution and fit 
projection. 



obtained. A second challenge is to determine 
background models from the upper B s mass side- 
band as well as from wrong sign B s — > Df tt + 
events. Additionally, we determine from simula- 
tion the mass distributions and momentum cor- 
rection of partially reconstructed modes. The 
main systematic uncertainties come from back- 
ground modeling. The measured lifetime is ct s — 
455.0 ± I2.2(stat.) ± 7A(syst.) (im. This is the 
best measurement in flavor specific decays and 
agrees very well with both CDF and DO measure- 
ments in B s — ► J /if) <j) decays discussed in section 
3.1. We note that the lifetime measured in fla- 
vor specific decays is not exactly the inverse of 
the decay width as measured in B s —> J/ip <fi de- 
cays. The actual relation between the measured 
quantity in flavor specific decays and the actual 
lifetime is: r(B° s ) FS = ±(1 + (ff) 2 )/(l - (f^) 2 ) 
which includes a small correction which depends 



8 



G. Giurgiu 



on the width difference AI\ It is also interest- 
ing to note that this measurement together with 
the lifetime measurements in B s — » J/ip <j> decays 
from both CDF [T5] and DO [2T] are larger than 
the 2007 world average [18] (1.41 ± 0.04 pa). The 
new world average B s lifetime will be much closer 
to the B° lifetime as predicted theoretically. 

4. Summary 

We have presented recent CDF results on B c 
and B s meson properties. The best B c mass mea- 
surement was performed by CDF in fully recon- 
structed B c — > J/-07T decays. Measurements of 
the B c lifetime in semileptonic decays agrees with 
the measurement from the DO experiment and has 
similar precision. These measurements provide 
useful information that can be used to improve 
theoretical models used to study heavy mesons. 
The recent measurements of the B s lifetime in 
both B s — ► J/ip<f> and flavor specific B s — ► D s it 
channels confirm the theoretical predictions that 
t s i=s Td- The measurements of the decay width 
difference Ar and the CP violation phase j3 8 are 
statistically limited but open the road for exciting 
results with full Tevatron data sets. 

I would like to thank the organizers of BEACH 
2008 conference for an enjoyable experience and 
my CDF colleagues and Fermilab staff for their 
dedication which makes these results possible. 
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